Completely Positive Maps with Memory 
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The prevailing description for dissipative quantum dynamics is given by the Lindblad form of a 
Markovian master equation, used under the assumption that memory effects are negligible. However, 
in certain physical situations, the master equation is essentially of a non-Markovian nature. This 
paper examines master equations that possess a memory kernel, leading to a replacement of white 
noise by colored noise. The conditions under which this leads to a completely positive, trace- 
preserving map are discussed for an exponential memory kernel. 
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The theory of open quantum systems deals with sys- 
tems that interact with an environment. This physi- 
cally realistic situation is of crucial importance in physics 
in general. The overall system-environment state is de- 
scribed by a density operator that evolves unitarily. Most 
often, one is interested in the system alone, which is de- 
scribed by a reduced density operator obtained by tracing 
over the environment degrees of freedom. The evolution 
of the system state is described by a completely positive 
map [3. 

The earliest prominent example of an open quantum 
system is the Weisskopf-Wigner theory of spontaneous 
emission describing the irreversible decay of an atom to 
its ground state. The standard approach is to make 
approximations that simplify the mathematics without 
changing the underlying physics. The vast majority of 
the theory is built on the Markov approximation. It is 
assumed that the correlation time between the system 
and environment is infinitely short so that memory effects 
can be neglected. This leads to the ubiquitous quantum 
Markovian master equation that describes the time evo- 
lution of the system density operator. This differential 
equation has a canonical form in terms of a Lindblad 
superoperator. The Lindblad superoperator generates a 
completely positive dynamical semigroup The 
solution to the differential equation for the system den- 
sity operator defines a completely positive linear map. 

Certainly, physical situations arise where correlations 
between the system and environment exist for a small 
yet finite period of time. Or the Lindblad superopera- 
tor may necessarily be of a more general form that de- 
pends on time. This complicates the theory, which can no 
longer be presented in the framework of quantum dynam- 
ical semigroups @ . To date, there lacks a solid theory of 
system-environment interactions in which memory effects 
are incorporated. In this Letter, we present a model that 
describes system-environment interactions with memory. 
Because this occurs in all branches of physics [|| it is 
of considerable importance to develop a general formal- 
ism that extends beyond the Markov approximation and 
resulting Lindblad evolution. 



We first review the canonical form of the generator 
that leads to the quantum Markovian master equation for 
the system density operator. The Gorini-Kossakowski- 
Sudarshan generator |jj has the form 



. n -1 

Lp = -i[H, P] + ^J2 ^{[FupFi] + [F iP , Fj]} 



(1) 



for p G A4 n [M. n denotes the set of n x n complex matri- 
ces), where H = H\Tt{H) = 0,Tr(i^) = 0,Tr(F i F J ) = 
Sij, and (cij) is a complex positive semidefinite matrix. 
The first term describes the unitary evolution while the 
second term defines the Lindblad superoperator C de- 
scribing the dissipative dynamics due to the interac- 
tion of the system with the environment. The solution 
p{t) — e Lt p(Q) to Eq. ( defines a linear operator 
<&t '■ P — y Pt that maps the system density operator at 
some initial time to the system density operator at some 
time in the future. The expression in Eq. Q has been 
shown to generate a completely positive dynamical semi- 
group that has the following properties: 

(i) || $ t p||i=||p||iVpe t>0, 
(a) $(®/„>ovn6Z + , (2) 

(Hi) lim uo $t = I, 

(iv) $ t $ s = $ t+s , t,s>0, 

where || • ||i denotes the trace norm in the space of linear 
operators on the Hilbert space H and V±~ (TL) is the cone 
of all positive semidefinite elements in the space . The 
first property states that the map is trace-preserving for 
all density operators for all times. The second property 
states that the map is not only positive, but completely 
positive, i.e., any extension to a larger space remains a 
positive map . The third property expresses the conti- 
nuity at the origin and hence for all time. Finally, is 
the semigroup property. This property arises when the 
environment or reservoir is delta-function correlated in 
time as is the case for white noise diffusive processes. 

The prevailing description for dissipative dynamics is 
given by p — Cp, where C is given by the second term 
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in Eq. QJ. Instead, we consider dissipative dynamics 
described by 

P = KC Pl (3) 

where K is an integral operator that depends on time of 
the form Kcj) = J* k(t - t')(j>(t')&t' . The kernel function 
k[t—t') is a well-behaved, continuous function that deter- 
mines the type of memory in the physical problem. The 
solution to the master equation can be found by taking 
the Laplace transform 

sp(s) - P (0) = K(s)Cp(s), (4) 

determining the poles, and inverting the equation in the 
standard way. The solution to Eq. © defines a linear 
map $ t : p — > p t that describes the evolution of a sys- 
tem coupled to an environment provided that <&t satisfies 
properties 2(z), (ii), and (Hi). Because the master equa- 
tion is no longer of the Lindblad form, the semigroup 
property is lost. However, this property is not necessary 
to describe a physically acceptable state evolution. All 
that is required is that the linear map be a completely 
positive, trace-preserving map. 

The evolution of quantum systems is described by uni- 
tary operators. Therefore, <& t should describe an evo- 
lution of the system that arises from an overall unitary 
evolution of the system and environment 

$ t (p) = Ti T {U(p®\ lo )( lo \)tf}, (5) 

where T denotes the environment degrees of freedom and 
7o is some initial state of the environment. One should 
imagine that the state p is prepared at some time t = 
and so is initially uncorrelated with the external system. 
The state and the environment evolve unitarily for some 
time and they become correlated. One may think of the 
environment as extracting information from the system 
as it will typically map pure states into mixed states. 
This noise process is described by a linear map involving 
only operators on the system of interest so that it has a 
Kraus decomposition 

<Mp) = £4m^ ( 6 ) 

k 

where the condition J^ fe = / ensures that unit trace 

is preserved for all time [jj . A map has a Kraus decompo- 
sition if and only if it is completely positive. Physically, 
the existence of a Kraus decomposition ensures that the 
system evolution is compatible with a unitary evolution 
on the system-environment Hilbert space. 

In solving Eq. @ it is advantageous, in practice, to 
find a damping basis |lOj | for the superoperator C that 
diagonalizes the master equation. Solving the eigen- 
value equation Cp = \p produces a complete, orthog- 
onal basis with which to expand the density operator 
at any time. This results in a set of eigenvalues {A^} 
and right and left eigenoperators, {Ri},{Li}, that sat- 
isfy the duality relation Tr {LiRj} = Sij. Once the ini- 
tial state is known p(0) = Ti{Lip(0)}Ri, the state 



of the system at any later time can be found through 
P(t) = E i ^{L i p(0)}A i (t)R l = EiTr^PCO)}^)^. 
The Ai(t) are general functions that determine the time 
evolution. The damping basis allows for the replacement 
of £ by the eigenvalues Xi so that Eq. (@J may be written 
as 

= E Tr{L J(O (0)}i?,. (7) 

i s — K [s)Ai 

We now present a concrete example of a system that 
has a master equation of the form of Eq. . Consider 
the time-dependent Hamiltonian 

3 

H(t) = fiJ2^i(t)(Ti, (8) 

i=l 

where Ti(t) are independent random variables and Oi are 
the Pauli operators. Each random variable obeys the 
statistics of a random telegraph signal, which is defined 
by Ti(t) — l) n< W. The random variable rii(t) has a 
Poisson distribution with a mean equal to t/2ri, while etj 
is an independent coin-flip random variable. By abuse of 
notation, the random variable ai takes the values ±a.;. 
The random telegraph signal is a wide sense stationary 
stochastic process |ll| with zero mean. This model is ap- 
plicable to any two-level quantum system that interacts 
with an environment possessing random telegraph signal 
noise. For example, this could describe a two-level atom 
subjected to a fluctuating laser field that has jump-type 
random phase noise . In the language of nuclear mag- 
netic resonance, this model describes a spin- 1/2 particle 
in the presence of three orthogonal magnetic fields. Each 
field has a constant magnitude at and inverts randomly 
in time with a distribution given by ni. The strength 
of the coupling of the system to the external influence is 
given by the parameters a^. The flipping or fluctuation 
rate is inversely given by t^. 

The equation of motion for the density operator is 
given by the von Neumann equation p = — j- [H, p] = 
—i J2k rfc(i)[o7c, p] that has the solution 

p(t)=p(0)-i[ J2r h (8)[* k ,p(8)]d8. (9) 

Ja k 

Upon substitution of Eq. © back into the von Neu- 
mann equation and performing a stochastic average, one 
obtains 

#*) =~ I H e^aliak, [a k , p(t')]]dt', (10) 
Jo k 

where the correlation functions of the random telegraph 

it-t' 

signal ( Tj(t)Tk(t')) = a\e T <= Sjk have been employed. 
Equation l|10|l is an incoherent combination of exact 
equations of motion for the mean value of the density 
matrix, i.e., the ensemble average over the random fluc- 
tuations for each component k = 1,2,3. After averaging 
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over the reservoir variables, we are left with a homoge- 
neous Volterra equation for the system density operator 
that has an exponential memory kernel. The state of the 
system at time t depends on its past history. The power 
spectrum of the environment is given by a Fourier trans- 
form of the exponential correlation function, which is an 
unnormalized Lorentzian with a maximum of 2a 2 r and 
a full width at half maximum equal to 1/(ttt). In the 
case of white noise, the delta-function correlation in time 
leads to a flat power spectrum for the environment with 
a strength given by a diffusion constant. The system is 
equally coupled to all frequencies of the external system. 
In the case of colored noise, the system prefers certain fre- 
quencies. Thus, a is the coupling strength of the system 
with the external system while r determines which fre- 
quencies the system prefers most. Increasing both a and 
1 /t corresponds to a more noisy environment. Therefore, 
the dimensionless product ar is the fundamental quantity 
that determines the amount of fluctuation. 

Having derived a master equation of the form of Eq. 
(PJ with an exponential kernel function, we now deter- 
mine the damping basis and solve the master equation. 
We assume that the fluctuation rates Tj are equal so that 
they obey the same Poisson statistics. This leads to a 
single kernel operator acting on the Lindblad superop- 
erator, rather than a linear superposition of such opera- 
tions: p — KiC\p + K2C2P + K3C3P. The damping basis 
for Eq. HI U| l is found to be the following set of eigen- 
values and eigenoperators: {Ao, Ai, A 2 , A3} = {0, —4(a 2 + 
oj}),-4(a? + a§),-4(a? + a2)} and {R , R u R 2 , R 3 } = 
{Lo, Li, L 2l L3} = {I, <7i, CT2, (73} which are self dual. Us- 
ing the damping basis, the Laplace transform of Eq. I|1(J|) 
becomes 



P(s) = £ 



s + 1/t 



s(s + 1/t) - A, 



Tr^O)}^. (11) 



This can be inverted to give the solution 



p{t) =^Tr{(7 i p(0)}A i (t)c7 i 



(12) 



In terms of the dimensionless time v 



tions Aj(f) 



cos(piis) 



- t/2r, the func- 
are damped har- 



monic oscillators having frequencies pi = \J (AniTi) 2 — 1 



with k 2 



a k for i 7^ j 7^ k. This differs from the 



Markovian case, where the functions are purely expo- 
nential functions in time with parameters defining the 
characteristic lifetimes. A power series expansion gives 
A(i/) = 1 - \{p 2 + l)v 2 + 0(iy) 3 , which shows that the 
linear term in v is missing. Thus, the standard white 
noise diffusion term vanishes. This is a general property 
of the memory kernel and a fundamental difference be- 
tween white noise and colored noise. 

The function A(z/) has two regimes - pure damping and 
damped oscillations. The fluctuation parameter, given 
by the product kt, determines the behavior of the so- 
lution. When < kt < 1/4 the solution is described 



by damping. The frequency p is imaginary with mag- 
nitude less than unity. When kt = 1/4 the function 
A(i/) = e~ u (l — v) is unity at the initial time and ap- 
proaches zero as time approaches infinity. In addition 
to pure damping, damped harmonic oscillations in the 
interval [—1, +1] exist in the regime kt > 1/4. 

The functions A»(i/) determine the evolution of each 
component of the Bloch vector. The bound |A(i/)| < 1 
ensures that the density operator evolves only to states 
on or inside the Bloch sphere so that <I>t always maps 
positive operators to positive operators. For Markovian 
master equations, the dissipation results in a contraction 
of each component, which is a consequence of the semi- 
group property. This property is absent for colored noise, 
so that the dissipation results in contractions with oscil- 
lations. The information exchange between the system 
and the environment leads to an exchange of entropy be- 
tween the two. The entropy of the average system state 
can oscillate in time with an overall decay. 

The solution to the master equation l|12(l defines a 
linear map <!>t : p — > pt on M2. This map has a 
Kraus decomposition $t{p) = J2k A k pAk with Kraus op- 
erators given by A 1 = y^iiyj en, A 2 = \/&(y) cr 2 , 

^3 = yiziy) °"3) an d A4 — \f^dy) I, provided the 
following linear combinations are nonnegative: 



4&(i/) = 1+Ai - A 2 - A 3 > 0, 

46(f) = 1 - Ai + A 2 - A 3 > 0, 

4&(i>) = 1-Ai-A 2 + A 3 >0, 

4£ 4 (f) = 1 + Ai + A 2 + A 3 > 0. 



(13) 



We now show which properties (2) hold. Clearly, this 
mapping is trace-preserving and property 2(i) is satisfied. 
Property 2(m) is satisfied because lim t | Ai(t) = 1. The 
system evolves continuously in time and the evolution 
is described by the identity map at the initial time. The 
semigroup property 2(iv) is lost as Aj(i)Aj(s) ^ Aj(t+s). 

The map is completely positive if <f>t®In > V n e Z + . 
It is sufficient to show that the composite operation on 
a maximally entangled state is positive [9|. For a linear 
map from Ai 2 to A4 2 -, we need only show that the com- 
posite map $t <g> I (|/3qo) (Poo I ) on M4 is positive semidef- 
inite, where |/3 o) = --^ (|00) + 1 11)) is the maximally 
entangled Bell state. Hence the map is completely posi- 
tive if and only if the eigenvalues of the composite map 
are nonnegative. The eigenvalues {£;} are exactly those 
given by Eqs. i|13|) . Therefore, property 2(ii) is satisfied 
for all a{T and every v if and only if 



inf X 3 (v) >0, 

(air) 



j = 1,2,3,4. 



(14) 



This follows from the fact that the existence of a Kraus 
decomposition and the composite operation on the max- 
imally entangled state are both necessary and sufficient 
to show complete positivity. 

Condition l|14|) is not satisfied for all values of the pa- 
rameters. A case where it is satisfied occurs when two 
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or more of the a, are zero. For example, if two of the 
dj are zero and one is nonzero then <I> t defines a de- 
polarizing channel with colored noise. Suppose a 3 = 
a,ai = a-2 = 0. In this case, Ai(z/) = A2(f) = A(z/) 
and A3(i/) = 1. This map has Kraus operators given by 
A a = V[l + AW]/2 / and v4 2 = - A(i/)]/2 <t 3 . 
This map is a completely positive, trace-preserving map 
for all values of the fluctuation parameter ar because 
|A(f)| < 1 for all times v. The Hamiltonian H = fiT 3 a 3 
implies that the z component of the spin- 1/2 particle is 
a constant of the motion so that the two states |±) a are 
fixed points and do not evolve. As a physical model, a 
constant magnitude magnetic field is applied in the z di- 
rection which inverts randomly in time. If the field is 
unobserved, the average trajectory for the density oper- 
ator of the spin- 1/2 system is dissipative. As time ap- 
proaches infinity, there is maximum uncertainty in the 
x and y components of the spin system. A steady state 
$>t(p)ss — \{p+o~ 3 po- 3 ) is reached; the entire Bloch sphere 
evolves to a line connecting the north and south poles. 

An example where condition l|14|l fails is given by 
ai = <22 = a and 03 = 0. Then each component 
of the Bloch vector evolves with frequencies given by 
Mi = M2 = \J (4ar) 2 — 1 and /i 3 = ^/32(ar) 2 — 1. There 
are no fixed points except the maximally mixed state. 
In this case, it is found numerically that ar < 0.8 leads 
to a completely positive map. This means that there is 
a limit on the environment power spectrum. The band- 
width may be large provided the coupling is made small 
and vice versa. 

We find that if two or more of the 0^ are nonzero then 
there are regimes for the fluctuation parameters where 
complete positivity is lost. Assume the frequencies Mi are 
equal. Then only the last eigenvalue in (fP3")l need be con- 
sidered. Setting the time equal to v = nj \i we find that if 
the frequencies are less than or equal to 7r/Zn(3) the map 
defined by Eq. (|12fl is completely positive for all time. 
Complete positivity is lost as the frequencies become 
large. The case of three equal frequencies sets an up- 



per bound. We have the following sufficient condition for 
(CHEJl to be satisfied: If [i* = max{Mi, M2j Ma} < 7r/Zn(3) 
then the map is completely positive for all time. 

We recover the Markovian master equation by letting 
t — » and a — > 00 in such a way that 2a 2 r becomes 
a constant. The random telegraph signal reduces to a 
Gaussian white noise in this singular limit and Eq. 11U|) 
becomes 

p(t) = - f S(t-t')J2Hr[<Jkd^p(t')]}dt'- (15) 

This master equation is local in time, leading to Aj(t) = 
exp(— jit) in Eq. H12fl with inverse lifetimes 7; — 4k 2 t. 
We point out that, even in the case of white noise, there 
are examples of maps that are positive but not completely 
positive fft. IT3I fl4|. For Markovian dynamics, the rela- 
tions qi3[l are satisfied if and only if ji < jj + jk for all 
permutations of the indices. 

By construction of our example, complete positivity 
holds in the white noise limit. We conclude that the 
loss of complete positivity, when the memory kernel is 
present, is a feature of the colored noise. The validity of 
Eq. (|3J has been questioned 0. This Letter shows 
conditions under which a memory kernel is valid and 
shows that such conditions can be achieved. We have 
presented an example of a completely positive, trace- 
preserving map that results from a system-environment 
coupling that is essentially non-Markovian. From this, it 
was shown that a memory kernel can be physically valid. 
The results reveal that interesting features arise from the 
colored noise of the environment, which are not present 
in the limit of white noise. White noise is an idealiza- 
tion of real noises and under certain conditions cannot 
be used. Thus, more work needs to be done in the study 
of more general noises that include memory effects. 
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